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Abstract 

We suggest an approach to constructing physical systems with dynamical character- 
istics of the complex analytic iterative maps. The idea follows from a simple notion that 
the complex quadratic map by a variable change may be transformed into a set of two 
identical real one-dimensional quadratic maps with a particular coupling. Hence, dynam- 
ical behavior of similar nature may occur in coupled dissipative nonlinear systems, which 
relate to the Feigenbaum universality class. To substantiate the feasibility of this concept, 
we consider an electronic system, which exhibits dynamical phenomena intrinsic to com- 
plex analytic maps. Experimental results are presented, providing the Mandelbrot set in 
the parameter plane of this physical system. 



PACS number (s): 05.45.-a, 05.45.Df, 05.45.Xt 



One of the rich and fascinating sub-disciphnes in nonhnear dynamics is the theory of iterative 
complex analytic mappings. A well-known example is the quadratic map 

Zn+l = ^ — (1) 

where the dynamical variable z and the parameter A are both complex. The set of parameter 
values defined by the condition that the iterations launched from the critical point z = do 
not diverge to infinity is the celebrated Mandelbrot set, perhaps, the most well-known example 
of a fractal [0, 0. Among other interesting objects in the field of complex analytic dynamics 
one can mention Julia sets - fractal basin boundaries of the attractor at infinity on the z- 
plane ||, period-tripling and other unusual bifurcation phenomena 0, ||, |^, Siegel discs - 
domains on 2;-plane, filled by closed invariant curves, which appear near fixed points at the 
moment of stability loss via irrational eigenvalues [0, ^, §. Although some nontrivial physical 
applications of complex maps are known (for problems like renormalization group approach in 
phase-transition theory and percolation theory 0, |I^), it would be interesting to find examples 
of nonlinear systems manifesting one or more of the above mentioned phenomena in actual 
dynamical behavior. 

In fact, complex analytic functions represent a very special and restricted class of maps. 
Indeed, the real and imaginary parts of f{z) must satisfy the Cauchy - Riemann equations. If 
this is not the case, the dynamics become drastically different |Tl|, |T^. This circumstance 
forces one to ask the principal question: Do the phenomena demonstrated by complex analytic 
maps have any concern to dynamical behavior of physical systems? Recently this problem was 
posed and discussed by Beck [|I4[. This author has considered motion of a particle in a double- 
well potential in a time-dependent magnetic field, and proved that under certain assumptions 
a complex analytic map may describe the dynamics of the particle. 

The aim of the present paper is to suggest some general approach to constructing physical 
systems with properties of dynamics specific to the complex analytic maps. 

Let us start from the complex quadratic map (|I]) and, first, separate the real and imaginary 
parts. Designating z = x + iy and A = A' + iX" we obtain 



Xn+i = X' - xl + yl, yn+i = X" - 2xnyn- 



(2) 



The variable and parameter change 

^ = x + I3y, 7] = X- /3y, (3) 

Ai = A' + f3X", A2 = A' - /3A", (4) 

where 7^ is an arbitrary constant, transforms the equations into a set of two coupled 
real logistic maps 

Cn+l = Ai - + e{^n - Vnf, Vn+1 = >^2 ' vl + ^(^n " Vnf ■ (5) 

Here e = {1 + /?^)/4/3^ plays the role of a coupling parameter. It is worth noting that coupling 
in these equations is of a very special kind: It may be interpreted as an equal simultaneous 
shift of control parameters in both maps at each step of iterations, which is proportional to the 
square of the dynamical variable difference. 

It is easy to see that for any selection of (3 the coupling parameter is positive, and exceeds 
1/4. Nevertheless, we may consider the coupled maps (|^) for arbitrary values of e. To motivate 
this, we turn to the generalized complex numbers [^, [T^, |T^. For pairs of real numbers {x,y) 



written as x + iy it is possible to define different consistent rules of arithmetical operations 
setting i"^ = a + ib, where a and b are some real constants. The case a = —1, 6 = 0, ori^ = —1, 
gives rise to usual complex numbers, a = l,6 = 0, 'i^ = l-to the so-called perplex numbers, 
and a = b = 0, i'^ = 0-to dual numbers. Any other selection of a and b appears to be 
isomorphic to one of these three cases, which are known as elliptic, hyperbolic and parabolic 
number systems, respectively. 

In the case of perplex numbers instead of Eqs. (|) we have Xn+i = A' — — y^, = 
A" - 2xnyn- Then the variable change (D, §) yields Eq. (g) with e = {(3^ - l)/4:f3^. The 
coupling parameter can take either positive or negative values satisfying e < 1/4. (In fact, by 
an appropriate variable change this case can be reduced to two uncoupled real maps.) For dual 
numbers x„+i = A' — x^, yn+i = A" — 2x„?/„. By means of (^), we arrive at the coupled 
maps (^) with 6 = 1/4- the same value of the coupling parameter for any (3. 

Diagrams in Fig. 1 are charts of the parameter plane (Ai,A2) for the coupled maps (|^) 
for several values of the coupling parameter. To obtain them we produce a large number of 
iterations starting from the origin ^ = 77 = at each pixel of some area on the parameter-plane 
and analyze the asymptotic behavior of the iterations. Divergence is marked by white, aperiodic 



behavior by black, and asymptotically periodic dynamics by gray; respective numbers designate 
the periods. At e = 0.5 we observe exactly the Mandelbrot set (rotated by 45° in comparison 
with its usual depiction). For 0.25 < e < +oo this set continues to exist, but as a distorted 
version of the standard picture. At e = 0.25 the set corresponding to the confined dynamics 
turns into a number of strips. For e < 0.25 it takes the form of rhombus-like structures (at 
e = this is a square). Similar metamorphoses from the Mandelbrot set to the rhombus-like 
object was noted earlier in Ref. for the quadratic map ([l|) considered for the generalized 
complex numbers. 

It is known that a single real logistic map represents a universality class, which is associated 
with the period-doubling bifurcation cascade and includes many dissipative nonlinear systems 
(forced nonlinear oscillators, Rossler and Lorenz equations, etc.) |]T^. It may be thought that 
taking two copies of a system, relating to this universality class, and maintaining the appropriate 
type of coupling one can arrange the dynamical behavior characteristic to the complex analytic 
maps. (It is supposed that the control parameters for the both subsystems allow an independent 
regulation.) 

To illustrate the feasibility of this approach we have elaborated a real electronic system, 
which reproduces the dynamics of the coupled logistic maps (^. 



We start with a specialized analog device suggested in Ref. [0 to study the dynamics of 
nonlinear systems represented by iterative mappings, in particular, by the logistic map. Our 
system (see Fig. 2) contains two pairs of the sample-hold cells (marked by dotted frames and 
figures 11, 12, 21, 22 on the diagram); one pair represents a single real quadratic map. Each 
of these cells consists of an analog switch and a capacitor. In the regime of picking a sample 
the capacitor is linked to the signal source via the switch, and accepts a charge up to a definite 
voltage. At some moment the switch breaks, and then the voltage on the capacitor remains 
constant - this is the regime of holding or storage. The voltage from the capacitor governs the 
operational amplifier of large input and low output resistance, so the charge of the capacitor 
remains practically constant. The used regime of the operational amplifier ensures equality 
of the output voltage to the input one. A set of two sample-hold cells is governed by two 
sequences of non-overlapping rectangular pulses: the switches Ku and K21 are opened while 
K12 and K22 are closed and vise versa. Multipliers A^^i and N2 ensure quadratic nonlinearity 



to obtain squared values of the voltages corresponding to rj^. Output of the operational 
amplifier D is the difference signal voltage. It is squared by the multiplier to produce 
the signal — r/n)^, then multiplied by e and added to the control voltages corresponding to 
the parameters Ai and A2 by means of the operational amplifiers Li and L2. The presence of 
three variable resistors , i?^, gives a possibility to regulate parameters Ai, A2, and e, 
respectively. 

Using an oscilloscope in the experiment we could distinguish either dynamics taking place 
in a restricted domain of the voltages, or the voltages jump to some distant values (analog 
of divergence in the mathematical model). For periodic regimes the periods could be easily 
determined (in units of the period of pulses, which control the switches) from the picture on 
the oscilloscope screen. 

Fig. 3 presents two examples of the topography of the plane of two parameters, namely, of 
voltages regulated by the variable resistors and i?^. The values of the coupling constant 
are chosen to be e = 0.1 and e = 0.5 respectively. 

The first diagram demonstrates a rhombus-like structure, similar to that in fig. Id. Some 
peculiarities of the experimental plot distinct from the computer model may be explained by 
inevitable deflections from perfect quadratic nonlinearities at the edges of the working interval 
of voltages. 

The second diagram for larger coupling shows a formation remarkably similar to the Man- 
delbrot set. Although in the experiment it was not possible to resolve extremely fine details of 
the structure, the location of all main leaves of the " cactus" are in excellent agreement with 
the computer generated picture. 

Experimental measurements at different values of coupling confirm that the set on the 
parameter plane {1/1,1/2) corresponding to the confined dynamics evolve in accordance with 
the results of numerical computations for the coupled maps (|^). 

It is worth emphasizing that what we deal with is a real physical object, effected by such 
factors as noise and technical fluctuations of voltages. The elements are not perfectly identical, 
the nonlinear function is not perfectly x^, and so forth. A substantial circumstance is that all 
these factors do not destroy the phenomena of complex analytic dynamics, which we observe. 

Perhaps, an electronic system could be built also to realize in a straightforward way the 



dynamics of real and imaginary parts of the complex variable governed by the map (|I]). However, 
our approach seems potentially more interesting because we state a direction for further search 
for dynamical systems manifesting behavior similar to that of the complex iterative maps. 

As explained, the quadratic map used as a basic element in our construction, must be re- 
garded as a representative of the wide universality class, which includes many realistic physical 
systems and their mathematical models (associated with the period-doubling bifurcation cas- 
cade) [1^. Hence, in the case of a properly arranged coupling between two period-doubling 
elements of any nature, one may expect the whole system to demonstrate the phenomena of 
complex analytic dynamics. 
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Figure 1: Charts of the parameter plane for the coupled maps (§) in dependence on the coupling 
parameter: (a) e = 0.5, (b) e = 0.3, (c) e = 0.25, (d) e = 0.1. Divergence is marked by white, 
aperiodic behavior by black, and asymptotically periodic dynamics by gray; periods are shown 
by respective numbers. 




Figure 2: Schematic representation of the electronic device corresponding to the coupled 
maps (^. Dashed frames show the sample-hold cells. Ku, K12, K21, K22 are the electronic 
switches controlled by sequences of the rectangular pulses. An, A12, A21, A22, D, E, Li, L2, 
Hi and H2 are the operational amplifiers, A'^i, N2, and are the multipliers. 




(*) 0*> 

Figure 3: Configuration of the set corresponding to dynamics in a restricted domain in the 
experiment with the electronic circuit. The chart represent the plane of voltages (f/i, U2) with 
f/i ~ 5Ai, U2 — 5A2 controlled by variable resistors Ri and i?^. The coupling parameter values 
are e = 0.1 (a) and e = 0.5 (b). The hatching on diagram (a) marks the domains of complex 
behavior 



